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ABSTRACT 

The problem of a laminate composite in presence of a 
crack located normal to the bond lines is considered. Stress 
analysis of the limiting case when the crack extends to the 
bond lines is carried out. Integral transform technique is 
used to formulate the problem in terms of a singular integral 
equation from which the power of stress singularity around 
the crack tip terminating at the interface is obtained. The 
singular integral equation is solved numerically and the 
effect of materfal property es on t" h e s t' r e s s — iTTt'efrs'i'ty - f'a"c't'o*r 
is calculated. 
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INTRODUCTION 


The problem of a laminate composite with a crack normal 
to the interfaces was studied by Hilton and Sih [1]. The com- 
posite geometry in [1] consisted of a cracked layer bonded 
between two half planes of different elastic properties. This 
is an idealization of a many-layered composite where one con- 
centrates on a single layer and approximates the effect of 
outer layers by prescribing some average elastic properties 
to the half-planes. 

The aim of this study is to reconsider the same problem 
and analyze the case when the crack has propagated through the 
layer and is just touching the interfaces. The limiting cases 
when the matrix is rigid or has zero modulus of rigidity corre- 
spond to that of a finite strip with constrained or free sides 
respectively. The method of solution uses the displacement 
expressions derived by Sneddon [2] where the problem of a 
finite strip with a central crack normal to the strip bound- 
aries is solved. First a crack problem will be formulated in 
terms of dual integral equations which will then be reduced 
to give a singular integral equation. This standard technique 
has been previously used by Erdogan and Gupta in various papers 
concerning crack problems [ 3 ] - [ 5 ] in composite structures. 

The method used by Sneddon [2] and Hilton and Sih [1] enables 
them to reduce the problem to a Fredholm equation of second 
kind. It also requires the crack tips to be away from the 
interface, thus limiting its applicability. 
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FORMULATION OF PROBLEM 


Consider a laminate composite in plane strain conditions 
consisting of a single layer of width 2h, shear modulus y 
and Poisson's ratio v bonded to two half planes having shear 
modulus and Poisson's ratio as y 2 and v , respectively (Figure 
1). A crack of length 2a (a < h) is located centrally along 
the x-axis. Only the symmetric problem will be solved in this 
paper. The skew-symmetric case can be treated in an identical 
manner. Subscripts or superscripts 1 and 2 will be used to 
refer to the layer and the half-planes, respectively. 


Considering symmetric normal tractions on the crack sur- 
face being the only loads [1], we have the following boundary 
and continuity conditions. 

Continuity conditions at x = h 


u x ( h ,y ) = u 2 ( h ,y ) ; v x ( h ,y ) = v 2 (h,y) 

= j v Qi,y) - ^yQi.y) 

Homogeneous conditions at y = 0 

a xy ( x 5 0 ) = 0 » I x I < h ; a£ y (x,0) = 0 , |x| > h 

v 2 ( x ,0 ) = 0 » | x | > h 


(1) 


( 2 ) 


Mixed boundary conditions at y = 0 

CT yy( x ’°) = “ ’ l x i < a 

(3) 

v 1 (x,0)=0 , a < | x | < h 
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The displacement field for the layer as derived by 
Sneddon [2] is a superposition of well-known transform solu- 
tions [3] for a body with x = 0 and y = 0 as planes of sym- 
metry and an upper half plane symmetrical about the y-axis. 
It may be expressed as 


u^x.y) = 


K -1 


f) w -j rs I 

- f J {fj T ['M 1 '!) \ — g ! (n)]s-inh(nx) 


+ xg x (n)cosh (nx) Jcosnydn 


0 » f, (?) K “ 1 t-,, 

1 / — F — (-4 ?y)e" ?y sin£xd£ 


* i z 


o °° I ■' i ■ 

v j (x ,y) = ~ f {- [f x (n ) + — 2 — (n) ]cosh(nx) 


K. +1 


where 


+ xg i (n)sinh(nx)}sinnydn 

? <f> , ( £) k +1 __ 

+ - / — j — + Zy)e ^ y cos£xd? 


x j = 3 - 4v^ 


The corresponding stress field is given by the equations 


°ix (x ' y) 


/ [f (n)cosh(nx) + nxg (n)sinh(nx)]cosnydn 

n 1 A 


9 00 r 

£ / (E) (I -5y) e"^ y cos£xd£ 
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ff L(x,y) o * r 

= \ / { ( n ) + 2g : (n)]cosh(nx) 

+ nxg i (n) sinhnxjcosnydn (5) 

O 00 £- 

- f J <j) 1 (^)(1+Cy)e"^ y cos?xd? 


J { [f (n) + g. (n)]si nh(nx) 

0 1 1 

nxg i (n)cosh(nx)}sinnydn 

“7 £y<!> U)e~^ y sin£xd£ 

71 0 1 

Condition a* (x,0) = 0 is identically satisfied by this repre- 
sentati on . 

Similarly, displacement and stress fields for the half 
plane can be expressed as 

u 2 ( x ,y ) = £ / {pj-[f 2 ( n ) + ~ 2 — 9 2 (n ) ] + xg 2 (n)}e _nx cosnydn 


xy 


(x,y) 


2P. 


2 

TT 


v 2 (x,y) = 
°xx (x ’ y) 

a vv (x>y) 

2M 2 

a xy (x > y) 

2^2 


( 6 ) 

9 00 1 K + 1 

f / {^Cf 2 (n) | — g 2 (n)] + xg 2 (n)}e _T1 sinnydg 

9 00 Y 

= - £ / [f (n) + nxg (n)]e _Tlx cosnydri 
n 0 

= | / [f 2 (n) + (nx-2)g 2 (n)]e" nx cosnydr) (7) 

0 

r \ CO 

= - £ I [f (n) + (nx-l)g (n)]e~ nx sinnydn 
n 0 z 2 


- 5 - 


Again this representation identically satisfies the homogene- 
ous boundary conditions (2). The unknowns f , g , f , g and 

1 ** 1 2 2 

4) i will be solved by using four continuity conditions (1) and 
the mixed boundary conditions (3). Continuity conditions (1) 
may be written as 

[f 2 (n) + ( + nh)g 2 (n)]e" nh = 

K l“1 

- [f^n) — g 1 (n) ]sinh (nh) - nhg x ( n ) cosh ( nh ) 

p n °° <f> , ( ? ) k -3 

- “ / p [2n 2 + 4 — (n 2 +£ 2 )]sin£hd£ 

77 o ( ^ 2 +n 2 ) 2 2 

[f 2 (n) - ( - nh)g 2 (n)]e" nl ’ - 

K +1 

[fj(n) + ~ 2 — g x (n) ]cosh(nh) + nhgj (n)sinh(nh) 

o2°° K +1 

+ pp / o C2C 2 + 4 — (n 2 +? 2 )]cos?hd? 

71 0 ?U 2 +n 2 ) 2 2 

( 8.1 

p [f 2 ( n ) + nhg 2 (n) ]e _rih = f x (n)cosh(gh) 

^ i 

~ 54>U) 

+ nhg (n)sinh(nh) + / w cos£hd£ 

1 77 o (^ 2 +n 2 ) 2 

- p- [f 2 (n) + (nh-1 )g 2 (n)]e" nh = [f^n) + 9 X (n)]sinh(nh) 

M i 

4n S 2 <}> U) 

+ nhg (n)cosh(nh) - - 11 / ~ sin£hd£ 

1 77 o U 2 + n 2 ) 2 


- 6 - 



Mixed boundary conditions (3) can be expressed as 


a; 


11 


(x ,0) 


2U 


\ / 4> 1 (?) (1 +?y)e"^ y cos?xd£ 

1 i m 

y->o + 


+ § / {[^(n) + 2g (n)]cosh(nx) 

0 

+ nxg i (n)sinh(nx)}dn = - , |x|<a 

3v K +1 ? CO 

9ir (x ’ 0) = " — ? / ^(Osin^xd? = 0 , a < | x | < h 


(9) 

( 10 ) 


Note that in (10) the boundary condition for displacement 
derivative instead of displacement* is used. This is done to 
have same dimensions of the equations (9) and (10). 

In order to reduce the dual integral equations (9) and 
(10) to a singular integral equation, we define a new unknown 
function G(x) as follows. 

3 v 

—1 — (-X-T 0 -) — =— G-(-x-)— ‘-j |-x-|-<-h (-li-)- 

From (10), it is clear that.G(x) = 0 for a<|xj<h. Inverting 
(10) and using (11), we obtain 

< +1 a 

- -i — <f>U) = / G(t)sin£tdt (12) 

c o 

From (8), eliminating f 2 and g , we find 


Displacement Vj is half the crack opening displacement. 
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f i (n) (x 2 + ^ )e nh + gi (n)[(K 2 + )nhe 7 


+ « 2 sinh(nh) + ~~ (e nh + < i e " n h ) ] = E x (n) 


^(nM ~ - l)e" nh + g x (n) [sinh(nh) - ( -p - l)nhe' 

M i p i 


+ ^7 (K i enh + e ' r ’ h)] " E >> 


where 


E x (n) = - ^ [D s (n) + Djn)] - ic 2 [D j (n> - D 2 („n 


E 2 (n) » ^ CD 3 (n) - D 4 (n)] + [D^n) + Djn)] 


4-, 2 00 , ( ? ) 

D (n) = ~~ / j cosEhdE 

1 11 » (? 2 +n 2 ) 2 


D 2 (n) 


7 1 / 5 - sinEhdE 

11 » (E 2+ n 2 ) 2 


D 3 (n) 


<3 00 <(> ( £ ) K - 3 

I 11 / — [2ri 2 + -4— (n 2 H 2 )]sin?hd 5 

22 ° (E 2 + n 2 ) 2 2 


o 2 00 4 > ( K ) K +1 

D (n) = / 9 [ 2 C 2 + -9 — (n 2 + C 2 )]cos?hd? 

' 12 0 E( E 2+ n 2 ) 2 2 


(13) must now be solved for f and g x and substituted in (9) 
to obtain an integral equation in 4> , or, by using (12), for 
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G(x). Combining (12) with the first integral of (9), a simple 
Cauchy kernel is obtained [ 3 ] - [ 4 ] . Second integral in (9) 
would involve integrals in (15) which, when combined with (12), 
give the following results. 

(k^+1 )D j ( n) = / G( t )o ( h- t ) e“ n ^ h-t ^ dt 

“ d 

(K +1 )D < n ) = J G(t) [n(h-t) - l]e-' l(h - t, dt 
-a 

( 16 ) 

(K i + l)D 3 (n) - - / G(t)[ n (h-t) + ]e- r| < h - t )dt 

(k +7)0 - (n) - / S(t)[ n (h-t) - ] e -i( h - t )dt + / G(t)dt 

-a c -a 

G ( t ) being an odd function in t makes the second term in D (n) 

4 

vanish. The details of above reduction may be found in Appen- 
dix A . 

Using (16), the singular integral equation from (9) can 
be expressed as 

/ dt + / G(t)K(t,x)dt 

-a I_x -a 

where 

K(t,x) = / k ( t ,x ,n)e" n ( h_t ^dn 
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k(t,x,n) = 


{1 - 

A 2 {1 


X 

2 


e 

[1 - 4 A n he _2nh - X X e" 4nh ] 

L 2 12 

( 3+2qh ) e 2n41 } - 2nxA i A 2 sinh(qx)e 
2n(h-t)}{cosh(rix)(3-2rih~X i e" 2ri * 1 ) 


[A cosh(nx) • 
-2nh + 

+ 2nxsinh(nx) } ] 
(18) 


K }J - K 11 
1 2 2 1 

U + K y 
2 2 1 


_ ^2 - 
2 + S y 2 


Kernel K( t ,x ) is a Fredholm's kernel for a<h as it is 
bounded for all values of t and x in (-a, a) is a<h. The singu- 
lar integral equation (17) is equivalent to the final Fredholm 
equation obtained in [1]. The present formulation of the prob- 
lem in terms of a singular integral equation has various advan- 
tages. A very simple numerical method is available to solve 
the equation [6]. However, the most important advantage of 
the procedure is that the case when a=h can be studied very 
effectively. The integral equation (17) must be solved sub- 
ject to the following s i ngl e - va 1 uednes s condition. 

a 

/ G(t)dt = 0 (19) 

-a 

The Case of a = h 

Referring to Figure 1, when the crack extends to the 
interfaces (case of a broken laminate), the Fredholm kernel 
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K( t ,x) in (17) is no longer bounded and contains point singu- 
larities at t = h and x = ±h. To extract these singularities, 
we need to study the infinite integral in the expression for. 
K(t,x). It turns out that the asymptotic value of the inte- 
grand k(t,x,n) as n + “ gives rise to these singularities. 

The part of the kernel contributing these singularities may 
be expressed as 

K (t,x) = 7 k co (t,x,n)e' n(h " t) dn (20) 

b o 

where, from (18), k ro is found to be 

k (t,x,ri) = e ^[A cosh(nx) + A (1 - 2ri(h-t)} • 

{(3 - 2nh)cosh(nx) + 2nxsinh(nx) } ] (21) 

Using the f ol 1 owi ngV resul t [7] 


J n m e -n(2h-t) f sinh(nx) '1 dn = d^_ J 0 -n( 2h- t )f s i nh (nx H dn 
o Lcosh(nx)J dt m 0 l_cosh(nx)J 


d m 

r 1 > 

rx V] 

(22) 


- ^2 2 ) 

dt 

(2h-t) - x ! 

' . 2 h - t~h 



the singular kernel K s (t,x) may be written as 


K. ( t ,x) = ( A + 3A J 


2 h - t 
(2h-t) 2 - x ' 


2 A 


[(2h-t) 2 - x 2 ]' 


[x 2 t - ( 2h-t) 2 (4h-3t) 


4(h-t){ (2h-t) 2 (2h 2 -ht-3x‘ ;: ) + x^ ( 6t/-3ht-x^ ) } j ( 23 j 


2 / > u 2 


{ ( 2h - t ) 2 - x 2 } 
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To analyze the behavior of the unknown function G(t) near the 
end points, we must consider the dominant part of the singular 
integral equation which can be expressed as 

h 2 

I( G(t)[ U + 2-i 4 M h+x > 2 U + 12X 2< h+x >a7 - \ + 3x 2 l r-T?h+7; 


+ jl^Jh-x) 2 - 12(h- x >37 - \ + 3x 2 ltTr kriT7 ]dt 


PfxfO+K,) 
— 4ir 


x | < h 


(24) 


G(t) is now assumed to have an integrable singularity at 
t = ±h which can be expressed as [8] 


G(t) 


H(t) H.(t)e 1riT 

(h 2 -t 2 ) Y (t-h) Y (t + h) Y ’ 


1 1 < h 


(25) 


where 0<Re (y) <1 and H(t) satisfies a Holder condition in the 
closed interval -h<t<h. Technique of determining y requires 
studying (24) near t = ±h [8, chapter 4] and has been treated 
"iTTllTtiri'l - !''?! - [-r]—["9"]— a-n-d— [-1-0-]- — e-ons-Tder— t-he— f-o-l-l-ow-i-ng— see- — 
tionally holomorphic function 


4>U) 



GUI 

t-z 


dt 


1 h f H( t)e 1Ti ^dt 
77 -h (t-h) Y (t+h) Y (t-z) 


(26) 


According to [8, chapter 4] 


<j>( z) - Hlh) L_ + * ( Z ) 

( 2 h) Y S1 "^ (z+h) Y (2h) Y s i nuy ( z - h) Y 0 


(27) 


where o ( z ) is bounded everywhere except at the end points ±h 
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where it has the following behavior 

H (±h) 

1 4> 0 ( z) | < — , Re(y ) < Re(y) (28) 

| z±h j o 

To consider (24) we need to reduce (27) for z = x, 2h+x and 
2h-x, obtaining 


4> ( x ) 


H ( -h ) cot-ny _ H(h) cotiry + .*# x % 
(2h) Y (h+x) Y ( 2h ) Y ( h- x ) Y 


x | < h 


cj) ( 2h+x ) 


H ( h ) _J + 

(2h) Y siniTy (h+x) Y 



h<2h+x<3h (29) 


4> ( 2 h - x ) = ^ + 4>*(x) , h<2h-x<3h 

(2h) Y sin-rry (h-x) Y 

Substituting (29) into (24) we obtain 


1 

( 2h ) Y s i n-rry ( h+x ) Y 


[H ( -h ) cosiTy 


- {4X 2 y(y+l)-12X 2 y-X i +3A 2 }] 


H ( h ) 

(2h) Y siniTy(h-x) Y 


[cosuy 


+ \ {4A 2 y(y+l )-12A 2 y-X i +3A 2 }] 


= P(x) 


(30) 


where P(x) contains all the bounded functions. 

Noting that for this solution of a symmetric problem G(t) is 
an odd function, hence, H(t) is odd, i.e., H(t) = - H(-t), 
multiplying (30) by (h+x) Y and substituting x = -h and then 
multiplying it by (h-x) Y and substituting x = h, we obtain the 
following characteristic equation for the determination of y: 

2cos7ry + 4A 2 (y-l)^ - (^ 1 + ^ 2 ) = 0 (31) 
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The equation is identical to the one obtained in [9]. The 
root of (31) satisfying 0<Re(y.)<l turns out to be a real con- 
stant for any material combination. If a<h in Figure 1, the 
only singular kernel will be the Cauchy kernel 1 / ( t - x ) . In 
that case, above analysis gives the characteristic equation as 

cot-fry = 0 , y = 1/2 (32) 

which is the well-known singularity at the crack tip in a 
homogeneous material away from any boundaries. Singular 
integral equation (17) now takes the following form. 

/ dt + / G(t)K(t,x)dt + / G( t) K f ( t ,x)dt 

-a t_x -a -a 


7Tp( x) ( 1 +K X ) 

4TT 


x|<a, a<h 


(33) 


K f (t,x) = / [k( t,x ,n) - k oo (t,x,n)]e" T1 ^ h ' t ^dn 
T 0 

K s (t,x) given by (23) is a Fredholm’s kernel for a<h and be- 
comes singular for a = h. k(t,x,n) and k oo (t,x,n) are given 
by equations (18) and (21) respectively and are such that 
Kf(t,x) is a Fredholm's kernel for a<h. 


Solution of the Integral Equation 


To solve (33), we first normalize the dimensions with 
respect to a, by introducing the following transformations 


t^ 

a 


y = 


_x 

a 


a = r 

o h 


G(t) = G(ax) = (j>(x) 


p(x) = p( ay ) = s(y) 


(34) 
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Hence (32) can be expressed as 


i 

/ 

- 1 


/ #(t) [ + a K s ( ax ,ay) + aK f ( ax ,ay) ]dx 

TTs(y) (1+« i ) 


4y x 


y| <1 - h * \ s 1 


( 35 ) 


As indicated in the previous section, we can express 4> ( x ) as 


<(>(x) 


A(t) 

0-t 2 ) y 


(36) 


where y is given by (32) if a Q <l and by (31) if a Q = l. (35) 
can be solved quite conveniently for a Q <l by using the method 
developed in [6], The method may be extended for a Q =l as 
described in [5, 9, 10], where the Gauss-Jacobi integration 
formulas are used. Thus, (35) may be approximated by 


N 

I 

j = l 


V (T j )[ 


x . -y . 
J 


a{K s (ax j . ,ay . ) 


K f (ax., 


ay,)}] 


its (y i ) (1 +k x ) 

4ir; 


(37) 


where 

p ( -Y , -Y) ( T ) = o , (j=l ,...,N) 

N J 

p(l -Yf 1 -y) (y ) = o 5 (i = l . ,N-1) 

N-l 1 

and A.'s are the corresponding weighting constants [5]. This 
J 

provides us with N-l linear algebraic equations for N unknowns 
( x • ) , j = l,...,N. The additional equation is obtained by 

J 

using the single valuedness condition (19): 


( 38 ) 


l A.if»(x.) = 0 
.1 = 1 J J 


\p( r.) are numerically evaluated from (37) and (38). We have 
J 

the derivative of the crack surface displacement (11) as 


3 v 


3 x 


1 ( x ,0 ) = G { x ) = 


a 2Y H(x) 

( a 2 -x 2 ) Y 


x <a , a< h 


(39) 


where 


H ( x) = i|>(f) 


Using the most common definition of the stress intensity factor 
as being the strength of the stress singularity ahead of the 
crack, we can write for the crack tip x=a as 


K = 1 im /2 (x-a) Y a (x ,0) 
x+a yy 


(40) 


It can be shown [9] that 


1im a vv( x ’°) = 

x-*a,(x>a) YY 


3 v . 


- 2y* 1 im (x ,0) 


x+a , ( x<a ) 


where 


(41) 


* } 

y 


z », 

t+t 


y y 
K 1 2 


[ 


a<h 


1 + 2 A ( 1 -y ) 1 "2A ( 1 -y ) 


S i nTry y x +k y , 


y +k y 
*2 2^1 


] , 


(42) 


a = h 


Hence the stress intensity factor K can be expressed as 


K = - 2y* (|) Y /2 H ( a ) = - 2/? y* (|) T MD (43) 
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NUMERICAL RESULTS AND DISCUSSION 


The numerical results are obtained for a uniform pressure 
distribution on the crack surface. Also, since the case of 
a<h has been treated previously* [1], results only for a = h 
are presented. Thus the input function is 

a yy ( x ’ o) = “ p( x ) = “ P 0 ’ M <h 

Figures 2 and 3 show the variation of the stress intensity 
factor with respect to the ratio of matrix-to-1 ayer shear mod- 
uli u / u for different Poisson's ratio combinations. These 
results are also presented in Table I. From this table and 
Figures 2 and 3, it is seen that the stress intensity factor 
ratio K/p Q h Y decreases with a decrease in • This value 

as y /y -* 0 ° corresponds to the semi -infinite strip fixed along 
x = ±h , where the power of stress singularity and the stress 
intensity factor depend only on Poisson's ratio of the layer. 

As y /y decreases, the power y of the stress singularity 

2 1 

increases, reaching a maximum when y ^ /y i ~ 0 ( i • e • > when the 
layer is rigid). In this case, a non-integrable singularity 
Y = 1 is obtained. Effect of Poisson's ratio on the stress 
intensity factor is also well depicted by these graphs and 
Table I. Keeping the same Poisson's ratio for the layer and 
increasing it for the matrix, increases the stress intensity 
factor. A reverse effect is noted if we increase the Poisson s 

Since a different method of solution was used in [1], spot 
check of a few points was made using the method presented 
here. The results were in good agreement with those in L 1 J * 
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ratio of the layer and keep the matrix unchanged. 


TABLE I 


Stress Intensity Factor K/p Q h Y 


^2 

v i = 

0.2 

v i = 

0.35 

Hi 

v 2 = 0.2 

v =0.35 
2 

v 2 = 0.2 

v =0.35 
2 

0.1 

0.249 

0.289 

mm 

mm 

0.5 

0.644 

0.836 

SB 

■US 

1 .0 

1 .000 

1 .343 

0.789 

1 .000 

2.0 

1 .611 

2.128 

1 .122 

1.408 

5.0 

2.881 

3.587 

1 .657 

1 . 945 

10.0 

4.084 

4.769 

2.042 

2.266 

30.0 

5.757 

6.183 

2.460 

2.572 

CO 

7.332 

7.332 

2.831 

2.831 


Figure 4 shows the crack surface displacement for two 
sets of material combinations. The left ordinate scale corre- 
sponds to Epoxy- A1 umi num combination. When the crack lies in 
e-p.o.x-v— Ta,v.e^. we ob serve much larger crack surface displacement 
in the center of the crack than when the aluminum layer con- 
tains the crack. Right ordinate scale corresponds to Aluminum- 
Steel combination in which the displacements are one order of 
magnitude smaller than Epoxy-Al umi num combination. 

Stress fields in both materials can be quite easily com- 
puted by using appropriate equations. Stresses in the neigh- 
borhood of the crack tip in such a case have been presented in 
[9]. Another question of interest is that if loads are 
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increased, would the crack tend to propagate into the material 
2, would it reflect back in 1 at some angle, or would debonding 
take place along the interface. A fracture criterion based on 
the "maximum stress" concept has been proposed in [9] and can 
very well be applied here. Validity of this criterion (or any 
other criterion) would have to be established by experimental 
s tud ies . 
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APPENDIX A 


Consider the integrals given in (15). Substituting from 
(12) and changing the order of integration we obtain 


4n 2 ? £cos£hsin£t 


(Kj+l ) (n ) = -2 / G(t)dt / 


o ( ? 2 +n 2 ) 


dc 


(tc+DD (n) = -2 / G(t)dt to / £to l n £ hs j n£ t „ 
1 2 « 11 » (5 2 +n 2 ) 2 


(k.+DD^h) = -2 J G(t)dt to J [ 


2n : 


o (5 2 +n 2 ) 


(fl-i) 


tc -3 , 

+ ] si n£hs i n£td£ 

2 ? 2 + r,2 


U.+DDjn) - -2 / G { t ) d t toi/ [ 


25 


* (5 2+ n 2 ) 2 


Ki+1 , 


2 5(5 2+ n 2 ) 


]cos£hsin£td£; 


— Us'fng— the— prope-rt-i-e-s— o-f— t-r-i-g-onome-t-r-i-e— ■fu-ne-t-i-on-s— and— the 

fact that G(t) is an odd function, the integrals can be reduced 
to the following form. 


(Kj+1 )D x ( n) 


a a n 2 00 

/ G(t)dt to- J 
a 0 


gsing(h-t) d £ 
(5 2+ n 2 ) 2 


(Kj+1 ) D 2 (n ) 


a n r, 00 

/ G ( t ) dt ^ j 
a 71 0 


£ 2 cos£( h-t) d £ 
(5 2+ n 2 ) 2 
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( A - 2 ) 


a n 00 

( |C 1 + 1 ) D 3 ( n ) = - / G(t)dt ^ / [ 

” d 


2iy 


0 (C 2+ n 2 )' 


K - 3 , 

— 1 ]cos5(h-t)d£ 

2 5 z +r) 2 


a 9 z «" 

K +1 > D .<"> - / G(t)dt J [ ^ 


2£ 


K +1 

1 t 


]sin£(h-t)d£ 


2 5 (S 2 +n 2 ) 

From the tables of Fourier Transform in [ 7 ] we have 


? sin^ydC = - 1 - (1 - e _ny ) 

o ?U 2 +n 2 ) 2n 2 


/ 


~2 sin^ydC = ye" ny 


o ( ? 2 +n 2 ) 

Therefore, by appropriate differentiation on both sides. 


obtain 


oo 


f 


: .1 . , co-s-g-y d-g— =---' -e— y - 
? 2 +n 2 2n 


7 — — o cos ^ yd ^ = 77 H - ny ) e" ny 

o U 2 +n 2 ) 2 4T1 

/ rrr^ cosKy ^ = k (1 + ^ )e_ny 

0 U 2 + n 2 ) 

Using (A- 3 ) - (A- 5 ) with y = h-t in equations (A- 
obtain the desired result in ( 16 ) as 


(A- 3 ) 

( A- 4 ) 
we 


(A- 5 ) 
, we 


-21 - 



(k +1)D (n) = J G(t)n(h-t)e" ri(h " t) dt 
1 -a 

( k +1 ) D (n) = J G(t) [n(h-t) - l]e~ n(h ~ t) dt 
1 2 -a 

a < - 1 - r» f h - 1 ^ 

(k +1 ) D (n) = - J G( t) [n ( h-t ) + — ]e nin ?; dt 
1 3 -a L 

a K, +1 

(k +1)D ( n ) = J G(t) Cn(h-t) - — — ]e dt + 

1 4 -a c 


( A-6 ) 


a 

/ G ( t) dt 
a 
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Figure 1. Geometry of the laminate composite. 
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Figure 3. Stress intensity factor vs. y 2 /y x - 
Layer Poisson's ratio = 0.35. 



Figure 4. Crack surface displacement for 
different material combinations 





